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Abstract 

J> ■ Neutrino telescopes detect the emission of Cherenkov light resulting from the 

tracks and showers of charged and neutral current neutrino interactions. These 
tracks or showers are reconstructed using a corresponding probability density func- 
tion (PDF) which depends on measured time and location of the detected photons. 
We call a composite event the mixed detector response due to the juxtaposition 
of more than one Cherenkov light source (track or shower). This paper presents 
the construction of a generic PDF corresponding to a composite hypothesis. This 
O composite PDF is therefore useful to reconstruct an arbitrary event topology and 

q , to favor or discard a given event topology hypothesis. 

^ : 1 Introduction 

> ■ 

^ ■ Neutrino telescopes focus on the reconstruction of single neutrino-induced muon tracks or 

single showers from neutral and charged current neutrino interactions. There are however 
events with a more complex topology due to a superposition of Cherenkov photons emitted 
from more than one track and/or shower. There are several reasons to develop the ability 
of reconstructing composite events (referred to as composite event reconstruction in the 
following), for instance: in a cubic km neutrino telescope, down-going uncorrelated muon 
(originating from distinct atmospheric showers) track events happen at a rate of a few 
per second at 2 km depth in seawater or ice; exotic channels may result in nearly parallel 
track events (production from neutrino-nucleon interactions of f pairs in supersymmetric 
extensions of the standard model p] or of muon pairs from micro black hole evaporation 
in models with extra-dimensions [2]); PeV r double bang events with significant conse- 
quences for particle physics (one expects 0.1-10 resolved events per year in IceCube [1]). 
Furthermore, the prolific number of shower-track events from muons undergoing catas- 
trophic energy losses or from neutrinos converting in the detector effective volume concur 
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with the composite hypothesis. The dedicated method of composite event reconstruction 
presented here may not be well adapted to the specific topology of muon bundle events: in 
which case, the fading light output along the bundle axis due to decreasing muon density 
is the important ingredient in the determination of the bundle properties (primary cosmic 
ray composition studies [5]). 

In section |2l single source reconstruction and likelihood definitions are briefly reviewed. 

Section [3] considers the case of a composite event topology and demonstrates the 
combination of source PDF's into a generic multisource PDF. In contrast to the single 
source PDF, the composite PDF explicitly depends on energy; the contributions from the 
sources to the detected light pattern is contingent upon the relative weights of each source. 
Thus, the reconstruction of the directional parameters of a composite event implies the 
reconstruction of the energy parameters. The composite PDF extends our current ability 
to reconstruct single topological events. 

The reconstruction of composite events is performed in Section H] with the presenta- 
tion of an extreme case, reconstructing simultaneously propagating parallel tracks in the 
IceCube array [6] . The reconstruction of an event is the minimization process of attaining 
the best likelihood given a topology hypothesis. We intend to demonstrate that the com- 
parison of likelihoods which correspond to double muon and single muon event hypothesis 
may be used to favor one topology over another. 

The normalized PDF describing the arrival time of a single photon originating fron 
source k at receiver j, given it has been detected, is written p^'^fa) = p\ ■ In the case 
of the detection of Nj photons at times {ti, ...,£#.}, the corresponding PDF is written 
p^\ti, tjv,)- It can be derived from the single photon PDF through the relation 
p( k ^(ti, t ) = Nj\ Y\f J p\ . The factor Nj\ reflects the indistinguishability of pho- 
tons (any of the Nj\ photon permutations leads to the same PDF and should therefore be 
summed up). 

The probability to detect n photons from source k at receiver j is written fn^ and 
will be assumed to obey Poisson law /«, = e _At / u ra /n!, where jj, = /j,( k >fi is the mean 
number of detected photons (oor hits or photo-electrons). 

The Pi s and /i's depend implicitly on the relative orientation and distance between 
the source and the receiver, the optical properties of the medium, the receiver efficiency, 
etc. The /i's additionally depend on the source intensity (related to the source energy). 

Consider now K sources and the detection of Nj photons at times {ti, tjv>} at 
receiver j. The corresponding normalized composite PDF 

p^ j} ^p^( tl ,...,t Nj ) 

will be constructed from the fn s and the p\ k s defined above. We also conveniently 
define the composite density p]^. = Nj\p^ used in the reconstruction formulas. 

2 Single source reconstruction 

For illustration purposes, we consider a source which can be parametrized by giving a 
vertex q, the incidence angles 9, <p and its energy E (e.g. infinite tracks or showers), that 
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is n d o { = 6. For a source with a different topology, e.g. an extended source or a finite 
track, the parameters should be defined accordingly and n± .{ acquire the corresponding 
value; this will not affect the results below. The reconstruction of the parameters con- 
sists in maximizing the log likelihood expression for an arbitrary hit time series for each 
receiver [7]: 



o 
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where j G {1, .., M} is the hit receiver index and nut = Yjj -^Yr The last term represents 
the information from receivers which have not detected any photons. This term can be 
dropped, when the calculation speed is a concern, but will consequently be accompanied 
by a degradation in the precision of the reconstructed parameters. We notice that eq. (pQ) 
can be split in two parts depending separately on energy E and directional ({q, 8, 0}) 
parameters. This leads to introduce the reduced log likelihood formulas: 

1 M 

In £ dir = V lnpgf (2) 

"hit - n d . .f r-f J 

and 

ln£ E = — — I > >./■■;■"■■ V h,ft n I (3) 




where n d o { take the corresponding values (respectively 5 and 1). Expressions (J2]) and (j3J) 
can be used to reconstruct the corresponding subset of parameters in a first approxima- 
tion, when calculation speed is a concern. However, best results are obtained by using 
expression (pQ), which more fully exploit the available information. 



3 Generic composite event topology reconstruction 

In the case of a composite event, recorded hits do not provide the information from 
which source they originate. So, in this section, we consider a case with an arbitrary 
number of sources and an arbitrary number of hits to each receiver, thereby describing 
the construction of the generic PDF p^'^ given the sets (k G {1, K}): 

{# j) }i<i<^.and{/f^} <n<7v r 
We also exhibit the log likelihood expression to use for reconstructing composite events. 



3.1 Two-source composite 

The simplest composite event consists of two sources and the detection of one photon at 
receiver j, the PDF is given by 



{2,j} 



f f (2J) , Aid) f (2j) J2,j) 

Jl Pi JO ' JO Jl Pi 



(4) 
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where we apply a superposition principle for the photon, which originates from source 1 
but not from source 2 with corresponding weights f9_ and /q 2 '"^, and vice- versa. The 
normalization is given by integrating over all hit timesj, meaning that the hit was actually 
measured: 



dhpl 2 ^ = 1 



For two photons, we have 



{2j} 



P2 



2 fM p™ f ™ + f™ (p™ p™ + p^ p™) 

, o f(l-i) f (2,i) J2,j) (2,j) 
i ^ JO J 2 Pi P2 



(5) 



where combinatoric factors appear to account for the indistinguishability of the photons. 
This is generalized for an arbitrary number of hits: 



iV, 



^{2,i} _ K {2,j} 
Pn, — n N 4 X 



\^1\ 3 m)\ f Nj _ m m. J m X -PnNj-mPirNj-m+l-PvN 



(6) 



m=0 



{^} 



where the sum over {vr^} is the sum over C(N,m) = N\/((N — m)\m\) different combi- 
nation terms. is extracted from the normalization condition 
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dti ...dt N ,p { ^f = 1 
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3.2 A'-source composite 

Expression is generalized for an arbitrary number of K sources in the composite: 



PNi 



K 



{K,j} 

N 4 



E 



A' 



AT, 



{ki,-,k Nj } 
Kki<K 



fc=l 



(7) 
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PNj'^ * s written as a sum of K Nj {source - hit} configurations instead of a permutation 
over hit indices as in eq. (EJ). Eq. (JTj) should be used also when K = 2 (instead of eq. ([6])) as 
it reduces the number of terms (thus reducing the computation time). The normalization 
is given by 

K 



N, 



-1 



K 



{K,j} 



e n/^„(E«! 

{ki,...,k N .} k=l 1=1 «=1 



(8) 



l<kj<K 

This formulation for a composite PDF is independent of the specific expressions for 
Pi and fn (different k may correspond to different source topologies). 



lr rhis is worth mentioning, as it has not been accounted for in the past, that experimentally, the 
recording time window for one event has a finite width. This fact should be adequately reflected in 
the hit probabilities fn^ and the basic PDF p^\t), which should be renormalized within this time 
window. 
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3.3 Likelihood formulation 



The calculation of the generalized log likelihood expression (TjQ) requires 1) using the com- 
posite PDF defined above and 2) setting appropriately for the increased number of free 
parameters na.o.f- That is, substituting the one source expressions for the hit probabilities 
and the photon arrival time PDF with the corresponding composite expressions: 



fV-J} _^ f{K,j} pnH P {i,j} _^ P {K,j} 

JO ~^ JO anQ Nj ~^ ' 



where f^ K ' j} = Yl k / (fcj) and 

P {K,jy _ (jAKJh-i jxjy 

Nj — K^Nj ) PNj 

are the M unnormalized PDF's which absorb the hit probabilities in their definition. The 

sets 

{fi k ' j) }ke{i..K}^d{p?' j) } ke{1 .. K} 

are determined in accordance with the underlying topology of the K sources defining the 
composite PDF. This reads: 

M 

ln£^ = 1 (]>>4f J} + E ln /o^' } ) 0) 



01^=0} 



The explicit dependence on energy with the inclusion of the /A terms in the def- 
inition of P^ K '^ forbids the splitting of expression ([9]) in two distinct components: the 
reconstruction of the directional parameters cannot be disentangled from the energy pa- 
rameter. This is an important difference between the single and the composite source 
likelihood formulations: in the former, while it is admissible to reconstruct the energy 
parameter separately from the directional parameters, c.f. eq. ([T]) (because the energy- 
dependent fn'^ in eq. (J7|) are absorbed in the normalization for K — 1), in the latter 
(K > 1), these factors are intertwined (the fn^ explicitly enter the definition of the 
Pjf K '^ preventing the splitting of \nC^ K ^). 

These notations have the advantage of unifying the log likelihood formulation: the 
eq. (P) is expressed exactly like eq. (pQ), the information regarding the actual number of 
sources in the composite is contained in }q K ^ and Pj^'^. The single source log likelihood 
formulation is therefore recovered by setting K = 1. 

For instance, in order to reconstruct two muon tracks, the expression (Q which depend 
on 12 free parameters should be maximized. In the case of two parallel muon tracks 
traveling simultaneously, this number is reduced to 9 ({q[, 0, (ft, E\, E 2 , d sep , rj} where 
d sep and rj are respectively the distance between the two muon tracks and the angle of 
rotation of one muon relative to the other). 



3.4 Inclusion of noise in the composite PDF 

When the considered topology hypothesis consist in more than one source, the individual 
noise of the receivers and the environmental noise levels must be included as a new source 
and cannot be included in the single source PDF. Thereby the number of sources is 
increased by one. 
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4 Use Case 



We discuss in this section a use case prompted by a study of the evaporation of a micro 
black hole into a muon pair [3] (see [2] for related micro black hole phenomenology studies). 
This use case also applies to other exotic channels (e.g. f pairs pQ). 

4.1 Motivations 

A fraction of the micro black holes produced in the neighborhood of the detector evaporate 
into PeV muon pairs, among other particles, which will travel 10's of km before reaching 
the detector with typical energies around 3 TeV. These events are characterized by nearly 
parallel and simultaneously traveling muons. Their incidence direction distribution is 
peaking horizontally. The rate of resolved events traversing the IceCube [6] detector 
will depend on the low energy gravity scale parameters (such as the number of extra- 
dimensions and the Planck scale) as well as on the detector resolution power (driven by 
the medium properties and the spacing of the detecting devices). The latter is the object 
of this study. The background rate of horizontal atmospheric neutrino-induced muons is 
approximately 100 km~ 2 yr~ x above 1 TeV, within one degree from the horizon (see e.g. 
Lipari [9] or Tig [I] for conventional atmospheric neutrino flux parametrizations). This 
use case is a challenge for the application of the composite reconstruction, because of the 
conjunction of simultaneity, parallelism and closeness of the muons. 

4.2 Event generation 

The AMANDA simulation software AMASIM [8], including the depth-dependent ice prop- 
erties, adapted for the IceCube geometry has been used to generate 1000 single muon track 
events and 1000 double muon track events. Subsequently the directional parameters were 
reconstructed by means of a first guess method (direct walk) [7J, providing the seed for 
the log likelihood reconstructions, which were performed under the single muon track and 
the double muon track topology hypotheses. 

Parallel muon pairs were generated 3 km away from the detector center with — 10 
TeV such that they mostly cross completely the detector active volume: -100 m < z < 
100 m, — 7r < < 7r, — 7r < rj < Ti, 5 m < d scp < 250 m, -200 m < xy < 200 m, where z is 
a vertical shift for the first muon from the detector center, xy is a horizontal shift for the 
detector center in the x?/-plane, is the incidence angle in the iq/-plane; d sep and rj were 
defined section (13.31) . 99.5% of the generated muon pairs did pass a majority trigger of 
24 fired receivers (photomultipliers) within 2.5 /is), each muon having a mean energy of 
3 TeV when reaching the point of closest approach from the detector center. 

The energy at the detector and incidence direction distributions of the generated events 
reproduce only the main features of the aforementioned detailed study. This simplifica- 
tion of the physical case nonetheless provides a useful framework with which to test the 
capabilities of the composite reconstruction. 
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4.3 Normalized PDF 



An explicit expression for the PDF [10J (also discussed in [7j) which can be used in shower 
and track reconstruction is 

K£,M) = f^) p ^~ le ~ pt (10) 

Given a detected photon at a distance source-receiver d, this expression allows to calculate 
the probability of a given time delay t (actual arrival time minus geometrical time) of the 
photon. £ = d/X, p = 1/r + c/A a b s , where the variables A, r, and A a b s characterize the 
medium. This is described in details [llj. Variant to this expression, accounting for a 
gaussian time response of the receivers [11] and for variable optical ice properties w.r.t. to 
the depth [12] were used for the reconstruction (it can be shown that using the average of 
the absorption and scattering coefficients between the photon emission location and the 
receiver is a valid approximation). 



4.4 Hit probability 

To calculate the relative probability of a photon from one track w.r.t. another, the hit 
probability must be calculated. First, we express the unnormalized version of eq. (fTOj) . 
that allows to calculate the probability for a photon to reach and to be detected at some 
distance d with some time delay t in an absorbing medium: 

1 1 ft \£ * (d+c ice t) 

P{£,t) = ^tI-) e"~ (11) 



r(0 t\rJ 

so that we can derive the light intensity at distance d by integrating ffTTT) : 

W(d)= / P{t,t)dt= e 1 Aab S (i + c . ceT / Aa6s )iA = w -d, (12) 
Jt=o L J 

The explicit depth dependence of p, P, W and subsequent derived quantities is omitted 
below. It must be noted that this depth dependence is not negligible for a medium with 
varying optical properties in general and for the medium where IceCube modules are 
deployed (South Pole ice [13]) in particular. 

The average number p(d, 9, E) of detected photons is proportional to W(d) (for an 
ideal receiver), with an absolute normalization factor a(9, E). The factor a accounts for 
the photomultiplier quantum efficiency and for the relative orientation to the track 9. The 
average number of photons is expressed as 



p{d, 9, E) = a(8, E)w d + b 



noise 



where 6 n oise is the detector /environment noise contribution (note that the split is approxi- 
mative and a may also depend on d). In the following we drop the ^-dependency (however 
this may not be a satisfying approximation for short distances between the source and 
the receiver, as Cherenkov photons have not been isotropized through scattering). 

Fig. [1] shows the distributions of the cumulated number of hits (left) and receivers 
(right) for all the single muon events in the sample w.r.t. the perpendicular distance to 
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the generated single muon track. Dividing the first distribution by the second results 
into the distribution of the average number of hits, shown Fig. [2J from which [i is now 
extracted^] by fitting /i(<f) = aw~ d + b. We obtain the proportionality factor a « 7, 
the noise contribution b ~ 0.008 and w ~ 1.045. When we calculate w from the usual 
medium PDF parameters [7], we get a close agreement (w = 1.044 using r = 557ns, A c g = 
25 m, A a bs = 98 m). This serves as an independent confirmation that the result above is 
correct. In general, a is a function of energy (or even a = a(6(d), E)). In this specific 
study, a is held constant, as the muons all belong to the same energy range. 



filtered hit distance all PMT distance 




Fig. 1: Filtered hit and receiver perpendicular distances to MC track. 

To reconstruct the energy (or more exactly to make a calorimetric measurement of 
the energy deposit in the detector), it would suffice to find a parametrization (using the 
Monte Carlo) for /j,(d) as a function of E as well 

The hit probabilities fn'^ can now be calculated from /i(<i), so that the normalized 
p\ 2 '^ can also be computed for the full likelihood reconstruction. 

4.5 Reconstruction and favored hypothesis 

Reconstruction has been performed using eq. (Q for single and double track hypotheses 
with fixed energy (i.e. fi = /i(cf)). A dedicated receiver response filter did consist in 
rejecting photomultipliers with more than four hits within 8 fis around the trigger time or 
with hit response inconsistent with a single photo-electron (within the same time window). 
The quality of this filter can be understood by looking at Fig. |2j At large distances (more 
than few 10's of m), the probability of having more than 4 hits is negligible. At short 
distances (less than few 10's of m), this probability is no longer negligible, but the fit no 
longer reproduces the data. This is mainly due to the fact that at short distances, hits 

Alternatively, can be extracted using the distribution of hit receivers instead of the distribution 
of the cumulated number of hits. 

3 The energy resolution log (E rec / Emc) for horizontal muons (using the MC described above) crossing 
the detector with energies between 1 TeV and 1 PeV is approximately 20% with this method. However, 
with the actual digital IceCube array, the energy can be reliably reconstructed up to much higher values. 
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Fig. 2: Resulting average number of hits per receiver w.r.t. the distance from Fig. ([T]). 
The hit probability can be extracted: p^{n > 1) = 1 — e~^, where jx{d) = aw~ d + b = 
aW(d) + b. The fitted parameters correspond to >50% and ~10% probability of having 
at least one signal hit at resp. 50 m and 100 m distances. 



can be the results of overlapping, and thus unresolved photo-electrons (the MC simulates 
the old AMANDA data acquisition system). This introduces a bias affecting the single 
and muon pair event reconstruction. However the bias is expected to be mild, given the 
muon energies. 

Some of the results are shown Fig. [3] and HJ 

— Fig. [31 the distribution of the difference between single and double muon log likeli- 
hood reconstruction hypotheses for single muon and muon pair events clearly demon- 
strates how an hypothesis can be favored against another. Notice the remarkable 
discriminating power between single muon and muon pair events, when d™^ > 50 
m. 

— Fig. |U the plots illustrate the correlation between reconstructed and generated d sep 
(left) and 77 (right). Though no first guess reconstruction was used to set initial 
conditions for d sep and rj, the correlation is nevertheless very strong (the correlation 
factors are respectively 97% and 87% for d sep and rj (d Bep > 50 m)). 

Evidently, seawater neutrino telescopes would improve the resolution to shorter d sep 
given the longer scattering length of the medium. 



5 Conclusion 

We showed how to construct a multisource PDF by combining single source PDF's and 
how to formulate the corresponding likelihood, thus enabling composite event reconstruc- 
tion. 
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Fig. 3: Distribution of the difference between single and double muon log likelihood 
reconstruction hypotheses In £ S i ng i e — In £ pa i r for single muon and muon pair events (with 
and without a selection on the reconstructed d sep ). 




Fig. 4: Left: Correlation plots of generated VS reconstructed d sep . Only events with 
ln£ pair < 12.5 are selected. Right: Correlation plots of generated VS reconstructed rj (for 
d > 50 m). Only events with ln£ pair < 12.5 are selected. 



Then using a IceCube simulation, we demonstrated that events were successfully re- 
constructed in the difficult case of close parallel muons. 

While the construction of this PDF is completely generic, it may not be necessarily 
effective at high energy (e.g. PeV double bang events) or very close distances, when the 
number of photo-electrons becomes large which in turn increases the number of combina- 
toric terms in the composite PDF (unless a compact approximation of ([7j) for large number 
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of photo-electrons can be found). However, it may reveal useful for the reconstruction 
of muons from distinct air showers and open experimental perspectives for various exotic 
searches. 
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